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Abstract 

Equations of motion are derived for the normal and the anomalous single-electron density matri- 
ces of a Fermi liquid using a time dependent finite temperature generalized coherent state (GCS) 
variational ansatz for the many-body density matrix. Self-consistent equations for the order pa- 
rameter A allow to investigate the interplay of Coulomb repulsion and pairing attraction in ho- 
mogeneous and inhomogeneous Fermi-liquids with spontaneously broken symmetry such as high 
temperature superconductors. The temperature of the Kosterlitz-Thouless transition to the two- 
dimensional superfluidity is calculated. 
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I. INTRODUCTION 



The BCS ansatz for the ground state wavefunction of superconductors predicts a gap 
in the spectrun^i originating from the pairing of two electrons with opposite momenta and 
spin projections. This pairing is induced by the effective attraction induced by electron- 
phonon interaction. The standard BCS description of superconductivity, as well as Eliash- 
berg's extension to incorporate Coulomb repulsion between electrons^i^ do not apply for 
strongly-correlated electrons such as in high-temperature superconductors {YBCO)^^ In 
two-dimensional {2D) superconductors (e.g., in cuprates)^, the phase transition to the su- 
perfluid state takes place at temperatures below the mean-field phase transition temperature 
for the appearance of the gapS*^, similar to the Kosterlitz-Thouless phase transition. The im- 
portant role played by electron-electron correlations in the local order parameter above the 
phase transition temperature, while the system is still in the normal phase, was discussed^. 

Superconductivity in strongly correlated electron systems is often described using effec- 
tive Hubbard Hamiltonians^. Taking electron-electron exchange and correlations into ac- 
count, the ground state energies and the collective excitation spectrum of superconductors 
have been calculated by Density Functional Theory (DFT) or its Time-Dependent exten- 
sion (TDDFT)S»2iOii In the Oliveira-Gross-Kohn (OGK) DFT equations^iSio^ the normal 
p(r, r') and the anomalous cr(r,r') density matrices satisfy the generalized Bogoliubov-de 
Gennes equation82J^ which include one external field Vs{r) coupled to the normal charge 
density n{r) = p(r, r) and a second field A(r, r') coupled to cr(r,r'). These fields contain 
exchange-correlation potentials fa;c(r) and A^dr, r'), respectively, obtained by requiring the 
charge density p(r) and the anomalous density matrix cr(r, r') of the interacting and nonin- 
teracting systems to be identical. TDDFT requires the same conditions to be satisfied for an 
externally driven system at all times^. The zero-order approximation in the standard DFT 
obtained by neglecting the exchange-correlation fields results in the Bogoliubov-de Gennes 
equations^ii^i which take into account the BCS pairing and classical Coulomb correlations at 
the Hartree level. 

The DFT equations for superconductors are usually formulated for the Kohn-Sham or- 
bitals in Hilbert space, and constitute a system of four self-consistent equations for Bogoli- 
ubov transformation coefficients u{r) and f (r), the normal density n{r), and the anomalous 
density matrix a{r, r')*. An extension of these equations to include a magnetic field was 
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proposed as wellii. 

In this article, we present an alternative derivation of the OGK equations for a Fermi 
liquid with spontaneously broken symmetry (FLSBS) based on a generalized coherent state 
(GCS) ansatz for the many electron density matrix. Coherent states were first used to 
describe anharmonic dynamical systems such as many body interacting fermions/bosonsi^ 
while preserving some of the useful properties of the original Glauber's coherent states 
for the harmonic oscillato r^'^i^^ . They encompass the Glauber coherent states as well as 
the squeezed statesi^^i^. The time- dependent Hartree-Fock-Bogoliubov (TDHFB) equations 
were derived for Boson systems using the GCS ansata^^*^. GCS are particularly suitable 
for variational dynamics by virtue the underlying Lie group algebr a^^i^^ . Using this ansatz, 
we derive equations of motion for expectation values of the normal and anomalous density 
matrices. We use the time-dependent variational principle which allows the description of 
the many-body system in terms of a small number of parameters and is formally closely 
related to classical Poisson bracket mechanics, i.e. to the variational equations of motion 
derived from the minimum action principle. A GCS representation for the BCS wave function 
has been used to analyze the coexistence of superconductivity and ferroelectricit y^^i^° . We 
obtain self-consistent equations of motion for the normal and anomalous density matrices^! 
and derive the spectrum of collective excitations, the density of the superfluid component 
at finite temperatures, and the temperature of the transition to the superfluid state for 
homogeneous and inhomogeneous superconductors. 

Reduced descriptions of many-body systems are naturally recast using density 
matrice o^^i^^ , and we found it useful to adopt the Liouville space density matrix 
represent ation22iSi of TDDF T^i^i^^ for the normal and anomalous density matrices. This 
requires solving only two equations of motion for the normal and anomalous density ma- 
trices p(r, r') and a{r, r') coupled to two artificial external fields. These fields contain an 
exchange-correlation contribution and guarantee the charge density and the anomalous den- 
sity matrix to be exact at all timea^. The ground state is the stationary solution of these 
equations of motion. This is in contrast to the four self-consistent equations for Bogoliubov 
coefficients n(r) and v{r) and the charge density n{r) = p(r, r) and anomalous density matrix 
cr(r, r')i2iS&. The density matrix is a two-point function compare to two one-point functions 
u and v^, nevertheless, the computational cost can be reduced. The reason is that the 
density matrices have non- vanishing elements only when |r — r'| is less than a coherence size. 
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which is typically very short. This allows to neglect many density matrix elements, making 
its size scaling linear rather then quadratic. The Liouville space representation provides a 
clear picture of the underlying coherence, since it is not possible to include the coherence 
size explicitly in the traditional Hilbert space computations. 

Using our GCS ansatz we further define expressions for first-order adiabatic (time- 
independent) exchange-correlation contribution to the order parameter, which include elec- 
tron exchange in the spectrum. This provides corrections to the charge density and the 
anomalous density matrix of Ref.— at each order in the perturbative series for the exchange 
correlation potential. The present approach is applicable for superconductors in general, but 
it is particularly relevant to strongly-correlated high-Tc superconductors, where the BCS and 
Eliashberg theories do not apply. 

II. THE COHERENT STATE FREE-ENERGY 

We start with the many-electron Hamiltonian H where the electron-electron interaction 
consists of both Coulomb repulsion V{v — r') = e^/|r — r'|, and pairing attraction W{y — v') 
between two electrons with opposite spins^ 

H = Y.j dr^lij) (-^V? - /i) ^.(r) 

+ iT^Jdr I rfr'V^t(r)V^t,(r')r(r-r')^.'(r')^.(r) 

vv' 

Here (r) and are the Fermi creation and annihilation field operators with the anti- 
commutation relations [■?/'(r), (r')]+ = (5(r — r') and [■?/'(r), ^/'(r')]^ = [■^/'''"(r), (r')]_|_ = 0, 
the indices v and v' denote the spin projections; m is the effective band mass of electron, [i 
is the chemical potential (Fermi energy), and e is the electron charge. 

We further expand the field operators in a single electron basis set 0ijy(r) 

V^t(r,t) = ^0*^(r)aL(t); 

i 

Mr,t) = ^<l),,{r)ai,{t), (2) 
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where and cii^, are the corresponding Fermi operators with the anti-commutation relations 
[aj,dj]+ = 5ij and [ai,aj]+ = [a|,a]]_|_ = 0; 0iy(r) are orthonormal atomic basis functions; 
and i runs over all basis electronic orbitals. 
Substituting Eq. (j2)) in Eq. ((T]) gives 



(3) 



i,j,k,l 



i,j,k,l 



Here tj,- is the single-electron matrix element 



2mh 



(4) 



Vijki is the Coulomb electron-electron repulsion 



ijkl 



j dridr2 0*( 



ri - r2 



r0fc(ri)0z(r2), 



(5) 



PVijfcz is an attraction responsible for the creation of electron Cooper pair. 

Eq. describes the interacting many-fermion system with Coulomb repulsion, and at- 
traction between two electrons with opposite projections of spin. In ordinary (BCS) su- 
perconductors this attraction originates from electron-phonon interaction^; in YBCO su- 
perconductors the short-range attraction results from the thermodynamically equilibrated 
phase ordering producing charge stripe order^, and it assumes in the following form 



W,jM = A / rfrirfr2 0*(ri)0*(r2) 



!)fclri )</)/( r2) 



(6) 



where n > 1 is a positive rational number; A and tq are constants, determined by system 
geometry. 

Our derivation is based on the following ansatz for the time-dependent many-electron 
density matrix. At zero temperature the system is in a pure state, and the density matrix 
is given by K{t) oc \ilj{t)){ilj{t)\, where the (unnormalized) many-electron wavefunction is 
assumed to be of the form 



\ip{t)) = exp 



exp 



ij vv' v,v'^v 



v,v'^v 



(7) 
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with l^o) being an arbitrary single Slater determinant^!. Eq. (|7j) is a generalization of the 
RPA and BCS wavef unctions: setting Ajj = it reduces to the Thouless representation of the 
single Slater determinants^, for hij = it reduces to the BCS ansatz for the superconductor—. 
At finite temperature T our ansatz reads 



K 



TreM-Ho/{kBT)) 



(8) 



where fc^ is a Boltzmann constant, and 



dt dr dr' 



Y: h{r, r', t)MrWAr') + h*{r, r', t)^t (r)^^,(r') 



+ Yl A(r,r',t)V^t(r)^t,(r')+ ^ A*(r, r', t)V^.(r)V^,,(r') 



(9) 



or using our basis set 



Y hij{t)dldj^' + Y hijit)^iua]^' + Y Aii(^)«tl«jV' 

uu' uv' u,v'^v 



+ Y ^ij{t)di,ydj^' 



(10) 



Here Tr denotes the trace in the many electron Fock space. A similar ansatz was recently 
used for Bose condensationi^. Eq. (jSI) constitutes a generalized coherent state (GCS)SSi2a (see 
Eq. ()A1|) ). Eq. ((Tj) is a limiting case of Eq. (jH)) obtained by a specific choice of parameters^; 
Eq. (jHJ thus holds at finite temperatures as well as at T = 0. 

The parameters h^j and Ajj will be determined variationally by minimizing the grand 
canonical free energy 



F(/i,T) = Tr{HK) - kBTTr{K\og{K)) =n-TS, 



(11) 



Here Ti. is the enthalpy, S is the entropy, and the chemical potential fi controls the average 
number of electrons through the following constraints!^ 



Tr (^KYald,}j =N. 



(12) 



Instead of using h and A as the variational parameters, we shall switch to the following 
variables: the normal density matrix 



(13) 



and the anomalous density matrices 



ai^j.^ = Tr{Kdii,dj_^), (14) 

(^iuj-u = Tr{Kdla]^^)., (15) 

The next step will be to express the free energy in terms of p and a. We start with the 
enthalpy. Since K is the exponent of a quadratic operator, we can use Wick's theoren>2i 
and express all averages of products of creation and annihilation operators with respect to 
K (denoted with a subscript 0) as averages of pairs of operators, which are generators of 
the closed algebra. In particular we have: 

{dla]dkam)o = -(alafc)o(a]am)o + (a]afc)o(aJam)o + (aJa])o(afcam)o. (16) 

Using this factorization we obtain for the enthalpy (see Appendix^, where for brevity we 
omit spin indices z/: 

H ^ Tr{HK) = W [hiP^J - P]^) + + A^a,,] . (17) 

Here /i is a matrix with elements 

hij = tij + —— '^{Vikij — Suu'Viikj)pkuiu- (18) 

k,i 

U is the volume and A is the order parameter matrix with elements 

Ajj = ^ ^ ^^ijmn^mn- (-^'^) 

mn 

We next require the expectation of the effective Hamiltonian (Eq- d)) to be the 
same as the expectation of H (Eq. Q). 

Tr{KH) =Tr{KHo). (20) 

It can be easily verified that the condition Eq. ()20|1 is met provided we set in Eq. (Iin|l 
hij = hij (Eq. (fTKjl ) and Aij = Aij (Eq. (fTn|) )^S. We next turn to computing the enthalpy. 
The effective quadratic Hamiltonian Hq given by Eq. (fTn|l can be alternatively recast in the 
form 

H = tr{QR), (21) 



where the matrix Q is defined as 

R is the generahzed single-particle density matrix 

r4 ' ' I. (23) 





and the symbol tr in Eq. ()21|) stands for the trace in the single electron space. 

Using the ansatz Eq. (jSJ, the many-electron system described by the Hamiltonian Eq. Q 
can be mapped onto the ideal system of non-interacting fermions with the quadratic Hamil- 
tonian Hq, determined by the matrix with the parameters defined by Eqs. (jl8p and (jl9|) . 
The entropy S of the system of fermions described by the quadratic Hamiltonian Eq. (pT|l 
is given hy^ 

S{p,a,fi,T) = -kBtr[f log f + {l - f)\og{l - /)], (24) 
where the matrix / is 

/ = ^ • (25) 

exp{Q/{kBT)) + l 

We are looking for the normal and anomalous density matrices p and a that minimize 
the free energy (Eq. (fTT|) ) together with Eqs. ()17|1 and ()24|1 assuming the GCS ansatz for the 
many-electron density matrix Eq. (jHj). This minimization yields the following equation^^ 

R = — i , (26) 

expiQ/{kBT)) + l 

which gives 

p = 1 - e[E]-^ tanh [E/{2kBT)) , (27) 
a = Ia[E]-^ tanh (^/(2A;bT)) , (28) 
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where the the order parameter matrix A is defined by its matrix elements Aj^ (Eq. (fT^ ). 

Eqs. f|?fj) , and constitute self-consistent equations for the equilibrium p , a and 
A. The E matrix is defined as 



E = VA2 + £^ (29) 



where matrix elements of e in the single electron basis set are 



kl 



Pkh 



and 



Substituting Eqs. (j2I|), ^ and ^ into Eq. ^ gives 



n{p, a, fi, T) = tr[2En + e - E + Aa - -Mp], 

2 



where the matrix elements of M 



kl 



Vikij — -Vikji 



Pkh 



and the n matrix is 

n = exp (E/2kBT) + l] . 
At zero temperature Eq. (j32|l gives the ground state energy. 



(30) 



(31) 



(32) 



(33) 



(34) 



III. EQUATIONS OF MOTION FOR GENERALIZED COHERENT STATES 

We are interested in the dynamics of the system coupled to two external fields: fs(r), 
which is coupled to the normal density p(r) = p(r,r), and Aext(r, r')) which is coupled to 
the anomalous density matrix cr(r, r')^*i2i. These fields account for exchange-correlation and 
provide a starting point for the TDDFT framework^. The total Hamiltonian then becomes 
Ht = H + Hext, where H is given by Eq. (jHl), and Hext represents the interaction with the 
fieldaSiia 

t 



Hext = E/ drVext{r,t)4jl{r)M^) - J J rfr' Ae.t(r, r', t)^t(r)^^,(r' 



(35) 



Our goal is to compute the dynamics of the system given by the total Hamiltonian Ht 
with the time-dependent external fields using the GCS ansatz for the many-electron density 
matrix Eq. (jH}. This can be accomplished by applying the closed equations of motion for the 



averages of GCS generators (see Eq. ()A9|) in Appendix These equations were obtained 
from the finite temperature time-dependent variational principle^ii^^. Since piuju, crluj-v 
and (J*^j^y are averages of GCS generators, we can immediately derive closed variational 
equations of motion for these quantities in real spac o ^>^. Substituting the parameters of 
the energy Ti (Eq. (fT7|)) and Hext (Eq- (ESI)) in the closed equations of motion for the averages 
of GCS generators Eq. ()A9|) . we get {h = 1) 

+ J dr2[V{r-r2)-V{r' -r^)] 

[p(r, r', t)p(r2, ra, t) - p(r, rg, t)p(r2, r', t)] 

- I dr2[Wir-r2)-Wir' -r^)] 

[a*{r, ra, t)cr(r2, r , t) - a{r, ra, t)a*{r2, r', t)] 

+ [Vextir,t) - Vext{T^',t)]p{r,r',t); (36) 

+ J dr2[V{r-r2)-V{r' -r2)] 

[a*{r, r', t)p(r2, rg, t) - cr*(r, rg, t)p(r2, r', t)] 

- I dr2[Wir-r2)-Wir' -r2)] 

[a*{r, r2, t)p(r2, r', t) - p(r, r2, t)a*(r2, r', t)] 

+ [Ae.t(r, r', t) - A,xt{r, r', t)](T*(r, r', t). (37) 

The order parameter A(r, r',t), which characterizes the excitation gap in the spectrum 
(Eq. (j2ni)), is defined as 

A(r, r', t)= J dr2Wir - r2)a(r2, r', t). (38) 

It is interesting to note that the stationary solution of Eqs. ()36|) and gives the normal 
p(r, r') and the anomalous density matrix cr(r, r'), which minimize the equilibrium free energy 
Eq. (fTT|) in the ground state (Eqs. (jTfjl and (j2HI))^'^- The calculation of the free energy in 
Sec. |n] is thus not necessary for the present derivation. The formalism of Appendix IX| 
allows us to proceed directly from the ansatz Eq. (jH)) to Eqs. and (jTTj) . The calculations 
of Sec. im provide a consistency check and connect our results with the more conventional 
derivations. 
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DFT is usually formulated in Hilbert space and involves the solution of four self-consistent 
equations for the Bogoliubov transformation coefficients u{r) and v{r), and the charge 
density n(r) = p(r,r), and the anomalous density matrix cr(r,r') (Eqs. (4.121), (4.129) 
and (4.130) in®). In the presence of the external time- dependent field, the ordinary TDDFT 
framework involves also the solution of a system of four self-consistent equations for Bogoli- 
ubov transformation coefficients u{r,t) and v{r,t), the local anomalous density cr(r, t) and 
the current density j(r, t) (Eqs. (20), (21) and (22) in^). Here, in contrast, we obtain the 
ground state free energy and density matrices p(r, r', t) and a{r, r', t) by the stationary so- 
lution of two equations Eqs. and (jHTjl . The solution of these equations gives the charge 
density n{r,t) = p(r,r,t) and the anomalous density matrix a{r,r',t). 

Eqs. (jSni) and ^7\i unify several widely used equations: if we neglect the pairing attraction 
W = and set cr(r, r',t) = o-*(r,r',t) = A(r, r',t) = the last integral in r.h.s. of Eq. (jHU)) 
vanishes, and Eq. (jHUj) reduces to the standard RPA equation^l. By neglecting the Coulomb 
repulsion V = 0,we obtain the BCS equationai (where in the r.h.s. of Eqs. (j36|) and pTjl the 
first integrals vanish). Neglecting the second term in the r.h.s. of Eq. (pTTj) gives Eliashberg's 
equationsi^, which incorporate into account Coulomb repulsion between electrons at the mean 
field level. 

Since our ansatzs for the many-electron density matrix Eq. (jHJ and wavefunction Eq. ((Tj) 
have the same number of variational parameters h and A, they yield the same equations of 
motion for the averages of the GCS generators (i.e. normal and anomalous single-electron 
density matrices). The only dependence on temperature and chemical potential is through 
the initial conditions (Eqs. (P7j) and (j2HI))- Eqs. and (pTTjl conserve the temperature and 
chemical potential at all times. 

In analogy with the RPA analysis^!, we can look for a solution for the density matrices 
in the form of the following equations for matrices 



Substituting Eq. (piJj) into Eqs. (^Bj) and ^7\. we obtain the spectrum of the collective 
excitations 



p = aX exp{—iLJt) + a*Y* exp{ii2i*t); 



a = aX exp{—iujt) + a*Y exp{iuj*t). 



(39) 




(40) 



where Skj and A^^ are given by Eqs. (^01) and (|^. respectively. 
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IV. APPLICATION TO THE KOSTERLITZ-THOULESS PHASE TRANSITION 



Recent studies of high temperature superconductors^i^li^ show a competition between 
two types of interaction between electrons: a pairing attraction, which makes the spectrum 
satisfy the Landau criterium of superfluidity by creating a gap in the excitation spectrun>i, 
and Coulomb repulsion which tends to eliminate the gap thereby destroying the supercon- 
ductivity. This competition leads to "stripe" high temperature superconductivity (HTS) 
in cuprates (e.g., Y Ba2Cu^0-j-s {YBCO)) above the low oxygen concentration threshold 



When A(r, r',t) 7^ the excitation spectrum satisfies the Landau criterium of su- 
perfluidity, and a three-dimensional (3D) system becomes superconducting. However, 
cuprates {YBCO) are two-dimensional {2D) structures'*^. The Kosterlitz-Thouless tran- 
sition temperature^ to the superfiuid state in a two-dimensional superconductive system 
is given by Tc = {Trh'^ns{Tc)) / {Aksmb),^''^' where ns{T^ is the temperature dependent su- 
perfiuid density of the superconductive system, and /c^ is Boltzmann constant. For tem- 
peratures close to the phase transition in the mean field approximation (Aoo(T'°) = 0), 
which satisfy T — T° ^ r°, the superfiuid density is ns{T) = (2(T° — T)n2Z))/ (^c );~ where 
^20 = PF/(2'/r^^) is the total two-dimensional {2D) density of electrons {pp is a Fermi 
radius). To find from the condition Aoo(7'°) = one needs to solve the self-consistent 
equations Eqs. (|^ - (^Hj) to obtain the temperature dependence of the order parameter A(T). 
The temperature of the phase transition can be estimated in the mean field approximation^^. 
For the gap spectrum of collective excitations Eq. (|4m) the mean field transition temperature 
T° is Aqo = 1.76/cbT°,— where Aqo is the order parameter at zero temperature. Combining 
these expressions, we obtain for the temperature of the Kosterlitz-Thouless transition, below 
which the superconductivity exists. 



Tc can be calculated using the order parameter obtained by solving Eq. (j2Hl)- 

Since both (Eq. dHJ) and the order parameter A (Eq. ()28|) ) decrease with the Coulomb 
electron-electron repulsion V (Eq. (jSJ) and increase with the electron-electron attraction W, 
Eq. (jni), Eqs. PT|l and (j28|l allow to study the interplay of Coulomb repulsion and pairing 
attraction between electrons in FLSBS'. Cuprates have two-dimensional {2D) structure^. 



{6 ~ 0.20). 




(41) 
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where the electron-electron Coulomb correlations may not be neglecte d^ The screened 
2D long-range Coulomb potential^ in the momentum space is V2d{p) = (27re^)/(p + K2d), 
where 2D Thomas- Fermi screening radius is density-independent = /?.^/(2mfee^). While 
the 3D short-range Coulomb potential in the momentum space V^dIp) = (47re^)/(p^ -|- 
^3d)) where 3D Thomas- Fermi screening radius^'' increases with density k^d ~ ng^^ and 
almost eliminates the Coulomb potential at large distances (r > k^^). Therefore, the 
Coulomb electron-electron correlations in Eq. ()28|) are much more important in the 2D 
order parameter compared to 3-D. The present theory thus describes the contribution of these 
correlations to the spectrum of collective excitations, the order parameter, superfluid density, 
the temperature of the Kosterlitz-Thouless transition and the DFT exchange-correlation 
potential. 

V. DISCUSSION 

Using the CCS ansatz (Eq. (jHl)) for the many-electron density matrix we have calculated 
the ground state free energy, the equations of motion for the normal and anomalous density 
matrices and the quasiparticle spectrum of superconductors. This results in a quadratic 
expression for the energy Eq. in the operators cij and a|, without using the Bogoliubov 
algebra of u — v transformations^. 

The advantage of the Liouville space representation of TDDFT for the normal and anoma- 
lous density matrices is that it only requires to solve two equations of motions for the normal 
and anomalous density matrices Eqs. (jHUj) - (jHTf) coupled to two artificial external fields fexi(r) 
and Aext(i", r'), which contain exchange-correlation^, instead of the four self-consistent equa- 
tions for Bogoliubov coefficients u and v and density matrices p and a. The normal p(r, r') 
and the anomalous cr(r, r') density matrices, which minimize the equilibrium free energy 
Eq. (jllj] ''^'^ in the ground state are simply given by the stationary solution of the Eqs. ()3(i|l - 
(jnZj)- And in order to get this ground state we actually don't need to derive the parameters 
of the effective quadratic Hamiltonian Eq. Q for the many-electron density matrix, as we 
did in SecHTl 

Our equations for the total energy, the spectrum of collective excitations and the gap are 
written in a general basis set and therefore apply to both homogeneous and non- homogeneous 
systems. In a homogeneous system we can use the plain waves basis, i.e. for the two- 
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dimensional system the eigenf unctions of a momentum p </'p(r) = f/~^/^ exp(— ipr); where 
U is the volume. Since we used a general basis set (not necessarily plane waves) our results, 
which contain non uniform normal p(r, r') and anomalous cr(r, r') density matrices, should 
be able to describe the short-coherence-length superconductors. For example, in YBCO, 
where the Inm coherence length is comparable to the lattice constant^. 

Finally, we comment on the connection of our results to the ground state energies and the 
collective spectrum of excitations in superconductors calculated using Density Functional 
Theory (DFT) and Time-Dependent Density Functional Theory (TDDFT)^'2ii2iii. Let us 
consider the following exchange-correlation potentials that depend on the density matrix 
(rather than merely on the charge density): 



Vxc[p]ir) = ^ ^ 



(42) 

n(r')=p(r',r') 



2 Sn{r) 

Substitution Eq. ()42j) in the DFT equations for the charge density and the anomalous den- 
sity matrijjSJii gives equations Eqs. (f?fjl and (f^H|) for the normal and anomalous density 
matrices. The GCS ansatz (Eq. ((Tj)) is thus equivalent to TDDFT provided we use the 
approximate exchange-correlation potential, Eq. ()42|) . To improve this functional, the adi- 
abatic (time-independent) exchange-correlation potentials can be obtained using the func- 
tional derivatives^ 

t;,.,([p,(T];r) = ^^^^^ ; A.,,i[p, a];r,r ) = ^, (43) 

where the exchange-correlation free energy F^c can be obtained using Feynmann diagram- 
matic perturbation theory for the self-energy of the Green function^LSk. The zero-order 
normal and anomalous density matrices derived using the present ansatz (Eq. (0)) are given 
by Eqs. and (jH7|) . respectively. The spectrum of the corresponding Green function is 
Eij given by Eq. ^ together with Eq. (IHUll . 

The first-order exchange-correlation contribution to the order parameter is identical to 
that of Ref.— for homogeneous superconductors provided we set Eij = Uj. The second term 
appearing in our expression for the energy Eij Eq. (p^ 



E 

kl 



Pkl 



comes from the first term aij{t)a\j^ajyi in the exponential in our ansatz Eq. (|7j), which 
represents electron-electron exchange. It corrects each order in the perturbative series of 
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Ref.— for the exchange correlation potentiaL This term is absent in the OKG equations 
which use as a reference the Bogohubov-deGennes approximation, which takes into account 
the BCS pairing and Coulomb correlations at the Hartree levelM and neglects exchange- 
correlation potentials'^. 
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APPENDIX A: GENERALIZED COHERENT STATE REPRESENTATION OF 
A HAMILTONIAN WITH ELECTRON-ELECTRON PAIRING 

Mathematically, a set of GCS is determined by a Lie group G, its irreducible unitary 
vector representation T with the space V and a reference state \Q) G V. The GCS are then 
states that have a form T{g)\Q) with g E G, where g is a set of parameter o^^i^^i^^ . 

Our ansatz for the many-electron density matrix (Eq. (jH))) can be expressed as 




(Al) 



where the set of numbers A, parameterizes the density matrix. The operator set {Tj} = 
{al^^aj±u, al^aj_^, ai^aj^u, 1} (/ is the identity operator) which forms the Lie group G is 
characterized by the commutation relations among the complete set of operators Tj necessary 
for describing the quantum dynamics of the system: 




k 



where G^j are known as the structure constants of the set {Tj}. Writing t/^ = diuCij-^, 
T^^^ = alj^a^j_j,, T-j'^ = a\j^^aj±i, + and, since dj satisfy Fermi anticommutation rules, 

we have the closed algebra of generators {T,} with respect to the following commutation 



rules 
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[fj), T;(+)] = 6,,6.^rT!^^ + 6,,,6.^,f^\ (A3) 

Eqs. ©, (EH), dH and (jUl) show that our ansatz for many-electron density matrix is 
generated by the operator set {Tj}, which forms the closed algebra of generators with respect 
to their binary commutation. The states described by our ansatz Eqs. (jHl) thus constitute 
generalized coherent states (GCS). 

The set of generators {Tj} corresponds to the Hamiltonian Eq. (^. The variational 
equations at zero temperature are derived as follows: given a Hamiltonian H, and time- 
dependent wave functions |f2(r)), we minimize the action: 

S[n{T)] = JdT[i {Q{T)\dQ{T)/dT) - {Q{t)\H\Q{t))] . (A4) 

By choosing \^{t)) to be a GCS, the resulting variational equations can be written in the 
Hamiltonian form for any set of coordinates Qj which parameterize 

^ = {n,n,} (A5) 

where {■ ■ ■} denote Poisson brackets and Ti is the classical Hamiltonian defined by: 

n{Q) = {Q\H\Q). (A6) 

The Poisson brackets clearly establishes the link between the variational equations and the 
classical dynamics. 

When the classical Hamiltonian is given by 

k 



n=l ii---i„ 



^ = E E C„(^n)---(T.J, (A7) 



the Poisson bracket assumes a very simple form provided the wave functions \Q) are parame- 
terized by the expectation values {Q\Tj\Q) of the operators Tj rather than by the parameters 
Qj. These expectation values then constitute a full set of parameters that uniquely specify 
the quantum state \fl). In particular, if Tj form a closed algebra Eq. ()A2|) . the Poisson 
brackets for (Tj) is given by: 

{{fj,{f^)} = zj2ctjn), (A8) 

A; 
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and the variational equations of motion for (Tm) take the closed form 




(A9) 



Using Eq. (|7j), it therefore suffices to derive the equations of motion for the expectation 
values {alttj) and (ojCij) to uniquely specify the dynamics of electrons. This may be done 
using the differential property of the Poisson brackets: 



When the expectation values are for generators of the set of GCS of some Lie group G, their 
Poisson brackets are given by the commutators of the underlying generators of the group. 
This direct correspondence between ordinary quantum mechanical commutators and the 
Poisson brackets greatly simplifies the calculation, since the variational procedure is then 
equivalent to the Heisenberg equations of motion. Eqs. (jHBj) and (jHTj) were obtained using 
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